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Abstract 

We examine a cosmological model with a dark energy density of the form Pde(0 = Px(0 + Pzif), where px is the 
component that accelerates the Hubble expansion at late times and pz(f) is an extra contribution proportional to H^{f). 
This form of pz{t) follows from the recent proposal that the contribution of zero-point fluctuations of quantum fields 
to the total energy density should be computed by subtracting the Minkowski-space result from that computed in the 
FRW space-time. We discuss theoretical arguments that support this subtraction. By definition, this eliminates the 
quartic divergence in the vacuum energy density responsible for the cosmological constant problem. We show that the 
remaining quadratic divergence can be reabsorbed into a redefinition of Newton's constant only under the assumption 
that ^(Oir^vIO) = 0, i.e. that the energy-momentum tensor of vacuum fluctuations is conserved in isolation. However 
in the presence of an ultra-light scalar field X with mx < Hq, as typical of some dark energy models, the gravity 
efl'ective action depends both on the gravitational field and on the X field. In this case general covariance only requires 
^^{T^y + {0\T^y\0)). If there is an exchange of energy between these two terms, there are potentially observable 
consequences. We construct an expUcit model with an interaction between px and pz and we show that the total dark 
energy density PdeCO = Px(t) + pziO always remains a finite fraction of the critical density at any time, providing 
a specific model of early dark energy. We discuss the implication of this result for the coincidence problem and we 
estimate the model parameters by means of a full UkeUhood analysis using current CMB, SNe la and BAG data. 
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1. Introduction 

Understanding the origin of dark energy is one of 
the most important challenges facing cosmology and 
theoretical physics (see e.g. [1-4]). One aspect of the 
problem is to understand what is the role of zero-point 
vacuum fluctuations in cosmology. In a Friedmann- 
Robertson- Walker (FRW) metric with Hubble parame- 
ter H{t) the bare vacuum energy density takes the form 

[Pbare(Ac)]FRW = [Pbare(Ac)]Mmk + O (//^(^A^) , (1) 

where [pbare(Ac)]Mink IS the bare vacuum energy den- 
sity in Minkowski space, whose leading divergence is 
(9(A|}), and we used for definiteness a momentum space 
cutoff Ac. In the usual treatment this A^ divergence is 
reabsorbed into a renormalization of the cosmological 
constant, giving rise to the cosmological constant prob- 
lem. The divergence oc H^A^ is instead absorbed into a 
renormalization of Newton's constant G [5, 6]. 

In this paper, expanding on results presented in [7], 
we reexamine the role of vacuum energies in cosmol- 



ogy. First, we wiU propose theoretical arguments sug- 
gesting that the correct way of computing the physical 
vacuum energy is to subtract the bare vacuum energy 
density of Minkowski space, [pbare(Ac)]Mink. from the 
FRW result given in eq. (1), before renormaUzing the 
result. By definition this subtraction eliminates the trou- 
blesome A^ divergence and, therefore, the cosmological 
constant problem. Then we turn our attention to the left 
over term H^A^ which now becomes the leading term 
in the vacuum energy. It is usually beUeved that this 
quadratic divergence can be reabsorbed into a renormal- 
ization of G. We show that this is correct only under the 
assumption that vacuum expectation value (VEV) of the 
energy-momentum tensor is conserved in isolation (an 
assumption that was implicit in the literature). However, 
general covariance of General Relativity (GR) only im- 
plies the conservation of the total energy-momentum 
tensor r^jy + (0|ruv|0), including both the classical term 
r^v and the semiclassical term (0|r^y|0). The separate 
conservation of (OITjuvIO) only takes place if we can de- 
fine an effective action which depends only on the grav- 
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itational field, by integrating out the matter degrees of 
freedom. This is possible only if the matter degrees of 
freedom are heavy with respect to the energy scale of the 
problem, and can then be integrated out. In a cosmolog- 
ical setting, this means that matter fields should satisfy 
m> Hq. If, in contrast, there is an ultra-light scalar field 
with m < Ho, as is typical for dark energy models such 
as quintessence, this field cannot be integrated out from 
the effective low energy action. We show that, as a re- 
sult, in general VT^^ = -V''<0|7;,vlO> ¥^ 0. In this case 
the effect of the quadratically divergent term in the vac- 
uum fluctuations cannot simply be absorbed into a re- 
normalization of Newton's constant G, and gives rise to 
interesting and potentially detectable cosmological ef- 
fects. We construct a specific coupled early dark energy 
model and test it against current observations. 

We use natural units where h = c - I, G - Mp^. If 
not specified otherwise, we work in a spatially flat FRW 
metric with signature ( — \-++), cosmic time t, scale fac- 
tor fl(f) and Hubble parameter //(f) - (da/dt)/a. Today, 
the Hubble parameter and the critical density take the 
values Hq and po = 3H^I{%nG), respectively. 

2. Subtraction of the flat-space vacuum energy 

In Minkowski space the divergence in the vacuum 
energy density is usually dealt with by normal order- 
ing the Hamiltonian, which gives by definition a van- 
ishing result for the physical vacuum energy density. 
However, it is useful to reahze that the problem can be 
treated more generally in the context of renormalization 
theory, which rather allows us to fix the renormalized 
vacuum energy density to any observed value. In the 
standard language of renormahzation, divergences in a 
generic A'-point Green's function are cured by adding 
the corresponding counterterms to the Lagrangian den- 
sity. The same procedure can be applied to vacuum 
energy, i.e. to the N = Green's function: one sim- 
ply adds a constant counterterm -pcount(Ac) to the La- 
grangian density. This corresponds to adding a term 
+Pcount(Ae) to the Hamiltonian density. Hence the re- 
normalized, physical vacuum energy density is given by 
Pren = Pbare(Ae) + Pcount(Ae). As always in renormaliza- 
tion theory, the counterterm Pcount is chosen so to cancel 
the divergences in pbare and leave us with the desired 
finite part that is fixed by comparison with the experi- 
ment. 

Using the language of renormalization theory is use- 
ful in this context because it makes clear that the cosmo- 
logical constant problem is not that quantum field theory 
(QFT) gives a wrong prediction for the cosmological 
constant (as it is sometimes incorrectly said). Strictly 



speaking QFT makes no prediction for the cosmologi- 
cal constant, just as it does not predict the electron mass 
nor the fine structure constant. Rather, it is a problem of 
naturalness, in the sense that the counterterm Pcount(Ae) 
must be fine-tuned to exceeding accuracy, in order to 
cancel the A^ divergence in pbare, leaving a physical vac- 
uum energy density that, if one identifies A^ with the 
Planck mass, is about 0(lO™) times smaller than A^. 

Posing the problem in terms of a cancellation be- 
tween Pbare(Ae) and Pcount(Ac) Can also give a first hint 
for a possible solution. First of all, one should ap- 
preciate that neither the bare vacuum energy pbaie(Ac) 
nor the counterterm Pcount(Ac) have a physical meaning 
and only their sum is an observable. Thus, this kind 
of cancellation is different from a fine-tuning between 
observable quantities. Indeed, the Casimir effect is a 
well-known example where a rather similar cancellation 
takes place. In that case the physical vacuum energy 
density of a quantum field in a finite volume is found by 
taking the difference between the bare vacuum energy 
density computed in this finite volume and the bare vac- 
uum energy density in an infinite volume. Regularizing 
with a cutoff both terms diverge as A^, but their dif- 
ference is finite and depends only on the physical size 
of the system. This might suggest that, similarly, to ob- 
tain the physical effect of the vacuum energy density 
in cosmology, one should compute the vacuum energy 
density in a FRW space-time and subtract from it the 
value computed in a reference geometry, which could 
be naturally taken as Minkowski space, leading to a sort 
of "cosmological Casimir effect". 

Before taking this analogy with the Casimir effect se- 
riously, one must however face the obvious objection 
that in special relativity the zero of the energy can be 
chosen arbitrarily, and only energy differences with re- 
spect to the ground state are relevant.^ In contrast, in 
GR we cannot chose the zero of the energy arbitrarily. 
One typically expects that "every form of energy grav- 
itates", so the contribution of Minkowski space cannot 
just be dropped. 

While it is certainly true that in GR the choice of 
the zero for the energy is not arbitrary, the point that 
we wish to make here is that what is the correct choice 
can be a non-trivial issue. As a first example, consider 
the definition of energy for asymptotically flat space- 
times. This is obtained from the Hamiltonian formu- 
lation of GR, which goes back to the classic paper by 



' Equivalently, one may observe that in the Casimir effect one actu- 
ally measures the force between the plates, i.e. not the energy density 
itself but only its derivative w.r.t. the size of the system L. Since the 
divergence is independent of L, it can simply be dropped. 
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Amowitt, Deser and Misner (ADM) [8, 9]. To prop- 
erly define the Hamiltonian of a given field configu- 
ration in GR one must work at first in a finite three- 
dimensional volume V. Then the Hamiltonian takes 
the form Hgr = //bulk + //boundary, where //bulk is given 
by an integral over the spatial volume V at fixed time, 
while //boundary is givcn by an integral over the two- 
dimensional boundary dV. When one evaluates //bulk 
on any classical solution of the equations of motion one 
finds a vanishing result (since //bulk is proportional to 
the constraint equations of GR), so the whole contribu- 
tion comes from the boundary term. On the other hand, 
//boundary diverges for any asymptotically flat metric g/^y 
(including flat space-time), when the boundary is finally 
taken to infinity. The solution proposed by ADM is to 
subtract from this boundary term the same term com- 
puted in Minkowski space rj^y. Accordingly, the energy 
E associated with a classical asymptotically flat metric 
g/jy is obtained by defining 

E = //grL?^,v] - ^^gr[v] • (2) 

This provides the standard definition of mass in GR, 
and reproduces the expected properties of asymptoti- 
cally flat space-times. For instance, when applied to the 
Schwarzschild space-time, it correctly gives the mass 
that appears in the Schwarzschild metric. This under- 
lines that our intuition that any form of energy gravitates 
according to GR is not entirely correct: eq. (2) tells us 
that the energy associated to Minkowski space does not 
gravitate. 

Similar subtractions also hold for non-asymptotically 
flat space-times, and can be performed either by sub- 
tracting the contribution of some reference space-time 
whose boundary has the same induced metric as the 
background under consideration [10-13], or even with- 
out introducing a reference background but just by 
adding some local counterterms to the boundary action, 
given by a coordinate-invariant functional of the intrin- 
sic boundary geometry [14, 15]. The latter prescription 
is particularly appealing for asymptotically AdS space- 
times. In fact, in the context of the AdS/CFT correspon- 
dence, this way of removing divergences in the gravita- 
tional action on the AdS side corresponds to the renor- 
malization of the UV divergences in the conformal QFT 
that lives on the boundary [14-18]. 

These examples show that, already in classical GR, 
the energy that actually acts as a source of gravity can 
be obtained from a Hamiltonian only after performing 
an appropriate subtraction. It is quite natural to assume 
that the same should hold at the quantum level, so in 
particular for zero-point fluctuations of quantum fields 
in curved space. To imderstand what is the appropriate 



subtraction for the FRW metric, we consider the Fried- 
mann equation that results from the Einstein equations 
sourced by Ty + (Oir^ |0), where Ty - diagi-p,p,p,p) 
is the ordinary contribution of matter, radiation, etc, and 

<0|r(^|0) = diag(-pvac, /5vac, Pvac, Pvnc) IS the Correspond- 
ing contribution of zero-point fluctuations: 

//'(0= ^(P+Pvac) . (3) 

We then require that Minkowski space. Hit) = 0, should 
be a solution in the limit p — > 0. This implies that 
all terms in [/Ovac]FRW that do not vanish for // — » 
must be subtracted. In other words, we must subtract 
the vacuum energy computed in Minkowski space. ^ In 
contrast, all terms proportional to or /Z'* are consis- 
tent with this requirement and thus generically allowed. 
This procedure eliminates the A^ term (as well as flat- 
space terms that appear for massive fields, such as m^A^ 
and m'' In A^, and are also much larger than the observed 
vacuum energy, for all known massive particles). Thus, 
this treatment of vacuum fluctuations solves the cosmo- 
logical constant problem by definition, at least in its 
"old" form, i.e. it explains why the observed vacuum 
energy is many orders of magnitude smaller than Mpj. 

We next turn to the quadratic divergence left over af- 
ter the subtraction. We will see that it could teach us 
something about the more recent form of the cosmo- 
logical constant problem, namely explaining why the 
observed dark energy (DE) density is just of the order 
of the critical density of the universe today, the coinci- 
dence problem. 

3. Non-interacting zero-point fluctuations and re- 
normalization of G 

Let us now discuss the fate of the quadratic diver- 
gence oc A^H^{t) in eq. (1). To understand what is 
the structure of the renormalized VEV of the energy- 
momentum tensor, consider the semiclassical Einstein 
equations for the renormalized quantities, 

^^V - \g>..R = »nG{T^y + {0\T^y\0)) , (4) 



^This conclusion also fits nicely with that of ref. [19] where the au- 
thors considered the QFT of a large number N of fields in Minkowski 
space, and found that the vacuum fluctuations collapse to black holes 

on scales smaller than O(N^^'^lfi). They conclude that Minkowski 
space would therefore be unstable to black hole formation unless ei- 
ther the length-scale where quantum gravity sets in, for a theory with 
fields, is of order A^'''*/pi, or the vacuum fluctuations in Minkowski 
space do not gravitate. 
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in which the vacuum expectation value of T^y is added 
as an additional source term. Together with the Bianchi 
identities the above equation implies 

v^(r^, + <oivio» = o. (5) 

This equation is therefore a consequence of the gen- 
eral covariance of the renormalized theory.^ If T^y 
and <0|r^y|0) are separately conserved we further have 
WT,y = V^<0| V|0> = 0. 

The stronger condition V(0|7),v|0) = can indeed be 
derived by using the effective action for gravity, which is 
obtained by treating the metric g^y as a classical back- 
ground and integrating out the matter degrees of free- 
dom (see e.g. refs. [21-23]). The VEV of the re- 
normalized energy-momentum tensor is then obtained 
by taking the functional derivative (2/ yf^)5l6g^y of 
the effective action. One can perform the calculation 
of the effective action using regularizations, such as 
dimensional regularization or point-sphtting, that pre- 
serve general covariance explicitly. The effective action 
is therefore explicitly generally covariant, and the VEV 
of the energy-momentum tensor derived from it is au- 
tomatically covariantly conserved. In this case, all di- 
vergences in <0|r^v|0) can be reabsorbed into generally 
covariant counterterms in the effective action. In par- 
ticular, the divergence oc H-^(t)A^ is cured by a coun- 
terterm proportional to the Einstein-Hilbert action. In 
fact, taking the variation of the Einstein-Hilbert term 
gives the Einstein tensor G^v, and Gqo, specialized to 
the FRW metric, is proportional to H^it). This shows 
that a VEV <0|roo|0) oc H^t) is obtained, in the ef- 
fective action language, by an additional term propor- 
tional to Jd'^xyf^R, and therefore is reabsorbed into 
a renormaUzation of Newton's constant, as it is well 
known [5, 6]. 

Such an effective action approach, however, assumes 
that we can integrate out all matter fields, i.e. that they 
are sufficiently massive with respect to the scale of in- 
terest. In our cosmological context this means that we 
are implicitly assuming that all fields have a mass m big- 
ger than the Hubble parameter H(t) at the time of inter- 
est, which is the quantity that fixes the relevant scale. 
Equivalently, we are assimiing that the wavelength 1 /m 
is smaller than the horizon size H~^{t). It is however 



'Observe that general covariance holds, for renormalized quan- 
tities, independently of the regularization scheme employed, as it 
should. If one uses a scheme that breaks general covariance at the 
level of regularization, such as a momentum space cutoff, then to re- 
cover it at the level of renormalized quantity one must also add non- 
covariant counterterms. See [20] for an extended discussion of this 
point, and comparison with the literature. 



interesting to consider the case in which in the spec- 
trum there is a scalar particle X with a mass nix < Hq. 
Such an ultra-light scalar field, with the mass protected 
against radiative correction by demanding that it is a 
pseudo Nambu-Goldstone boson, provides in fact a typ- 
ical realization of quintessence [24]. In this case we 
cannot integrate out this field, and the low-energy effec- 
tive action necessarily depends both on the metric and 
on this scalar field. Then the above derivation giving 
V<0|r^vlO> = no longer goes through. By taking the 
variation of this action with respect to the metric we 
get the total energy-momentum tensor, including both 
(0|r^v|0) and the energy-momentum tensor of this 
scalar field, and general covariance now only implies 

v'^«0|r^,|0) + T^^y) = 0. 

Whether these two terms are separately conserved 
is now a dynamical question, and depends on whether 
there is an interaction among them. One can imag- 
ine mechanisms by which vacuum fluctuations can ex- 
change energy with other forms of matter. Typical 
examples are the amplification of vacuum fluctuations 
[25, 26], or the change in a large-scale scalar field due 
to the continuous horizon-crossing of small-scale quan- 
tum fluctuations of the same scalar field, which is also 
at file basis of stochastic inflation [27]. If ^(Oir^vIO) = 
-V^T^fy + 0, it is no longer possible to reabsorbe the ef- 
fect of (0|r^v|0) into a renormalization of the Einstein- 
Hilbert term, nor of any other generally covariant local 
operator in the effective action. In fact, taking the func- 
tional derivative (2/ of a generally covari- 
ant term, we necessarily obtain a covariantly conserved 
tensor, so we can never obtain a quantity (0|r^y|0) that 
satisfies V''<0|7;,y|0) + 0. 

4. Cosmology with zero-point fluctuations 

We now explore the cosmological consequences of 
the hypothesis that (0|r^v|0) and T^v are conserved in 
conjunction, as in eq. (5), but not separately. In this 
case, as shown above, the term A^//^ cannot be reab- 
sorbed into a renormalization of G, and we rather ex- 
pect that it will give a genuine contribution to the to- 
tal energy density pzii) - 0{H^{t)M^), where M is 
the UV scale where new physics sets in (so that M 
could be typically given by the Planck mass Mpi, or by 
the string mass). Recalling that the critical density is 
Pcit) - 3H^{t)MlJ(8n) we see that, for M of order Mpi, 
pz(t) is of order of the critical density pdt) at any time 
t. Thus we write 

pzit) = ^zPcit) = QzPoHHt)/Hl , (6) 
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where po is the present value of the critical density. The 
value of Q.Z is fixed by the renormalization procedure to 
the observed value, so we keep it as a free parameter 
The same is true for the equation of state (EOS) param- 
eter wz defined by pz - wzPz, that can in principle be a 
function of time (see [20] for an extended discussion of 
this point). 

The late time acceleration of the Universe cannot 
be explained only by a DE density that scales like 
H^{t) [7, 28, 29], basically because observations tell us 
that the total DE density is at least approximately con- 
stant in the recent cosmological epoch. Therefore we 
assume pz(f) only to provide a part of the total DE den- 
sity that we write as pDE(f) - px{t)+Pz(t)- Here px(t) is 
a second dynamical DE component that is the one dom- 
inating the energy budget at the current epoch. In this 
approach the physical origin of px and pz can a priori be 
completely different. For instance, px could be due to 
a scalar field, as in quintessence models. In particular, 
we will take px to be due to an ultra-light scalar field 
with m < Hq as discussed in sect. 3. We define the EOS 
parameter of the Z-component by wx = Pxipx- Since 
Px and Pz have different origins, wx and wz can in prin- 
ciple be different and we will see that only for wx + wz 
a tracking mechanism for DE emerges. 

Mechanisms such as the amplification of vacuum 
fluctuations can produce an interaction between pz and 
scalar fields, but do not lead to interactions with pho- 
tons nor massless fermions, that are not amplified in 
the FRW space-time because of conformal invariance. 
As we saw in sect. 3, massive particles with m > Hq 
can be integrated out in the effective action and there- 
fore cannot contribute to the violation of the condition 
V(0|7'^v|0) = 0. Thus, in our context it is natural to 
consider an interaction between pz and px, while non- 
relativistic matter and radiation are conserved in isola- 
tion, so they scale in the standard way, pM ~ cT^ and 
Pr ~ ot'^ . In contrast, pz and px satisfy the coupled 
conservation equation 

Pz + Pa- + 3(1 -H wz)Hpz + wxWpx = . (7) 

In terms of the total DE density, PdeCO - Px(t) + Pzit), 
it reads 

Pde + 3(1-1- wx)Hpi)E - Mwx - wz)Hpz ■ (8) 

We insert pz(0 - ^zPoH^ / on the r.h.s. and use 
the Friedmann equation, that now reads H^(t)/Hl = 
Qfia"^ -I- Qjvffl + Pde(0/Po- With the definitions 

wo = Wx + e, e = Q.z{wz - Wx) , (9) 




Figure 1 : The energy fraction PDE{'^)/piot('^), for wq = - 1 and e = 0. 1 
(solid line) and e = -0.1 (dashed). 




Figure 2: The equation of state parameter woEio) for the same models 
as in fig. 1. 



and using d/dt = aH{d/da), we obtain 
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Note that the background evolution of the total DE den- 
sity is fully determined by the two parameters wq and 
e which can in principle be functions of time. We as- 
sume these functions to be constant for the scope of this 
work, i.e. we assume wz and wx to be constant. A more 
general analysis will be presented in [20]. Then the evo- 
lution can be solved analytically. 
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where C is the integration constant that is fixed by the 
condition pDE(a - l)/po - ^de = 1 - (^^r + ^m)- Ob- 
serve that the model reduces to wCDM when wx = wz 
and the deviation from wCDM only depends on e, so we 
are dealing with a one-parameter extension of wCDM, 
that we will caU wZCDM. 
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A very interesting feature of (11) is that poE always 
scales as the dominant energy component, for e 9^ 0. In 
fig. 1 we show the ratio pDE(fl)/Ptot(fl), where ptot(fl) - 
PDE(fl) + poQm«^ + Po^^R« Deep into the RD phase, 
as well as in the MD phase, pDE(fl)/Ptot(a) = (5(e), while 
today it becomes 0(1). Compared to ACDM, where the 
ratio Pa/ Plot is of order one today, but goes to zero as 
during MD and as a** during RD, the coincidence prob- 
lem is sensibly alleviated. 

We cannot claim that wZCDM entirely solves the co- 
incidence problem; in this model, in fact, the transition 
between the regime where Pde/Pioi = 0{e) to the regime 
where Pde/Pioi -0(1) takes place at the present epoch 
simply because we have fixed the integration constant C 
in (11) by the requirement that pDE(fl = l)/po = ^^de- 
Nevertheless the coincidence problem is certainly alle- 
viated, compared to ACDM, where poE is parametri- 
cally different from ptot, and the ratio Pde/Pioi evolves 
from (9(10 '^°) at a Planck time to ~ 0.7 today. 

It is also interesting to note that our model provides 
a different theoretical justification for parameterizations 
of early DE models that have been proposed in the liter- 
ature. To make the relation explicit, it is useful to define 
Wde = Pde/Pde and derive it from the total DE conti- 
nuity equation, pde + 3//(pde + Pde) = 0. Inserting the 
expUcit solution for pDE(f) found in (1 1) we get 

H'DE(f) = WO + e[pM(t) + PR(t)]/pDE(t). (12) 

This function is shown in fig. 2, setting wo = -1 for 
definiteness. For e > it evolves smoothly from a value 
wde -1/3 during RD, to wue - during MD and fi- 
nally goes asymptotically to wue - wq. This EOS is 
quite similar to that obtained in a commonly used pa- 
rameterization of early dark energy [30, 31]. Note that, 
for e < 0, wde goes through infinity, as a consequence 
of the fact that pde goes through zero with poE ^ 0, 
see fig. 1, but the pressure wdePde stays finite, and the 
background evolution is regular. 

Our model has some similarities, as well as important 
differences, with other DE models studied in the litera- 
ture. In particular, in [32-35] a model was proposed 
where the DE density has the form pA(f) = n() + n\H^{t), 
inspired by the idea that the cosmological constant 
could run under renormalization group. In this case, 
however, the two components no and n\H^{t) neces- 
sarily have the same EOS parameter. More closely re- 
lated is the AXCDM model proposed in [36], in which 
DE has two components, an energy density p^ - hq + 
riiH^it) associated to arunning of the cosmological con- 
stant interacting with an unspecified dynamical "cos- 
mon" field, although in our case the interaction is rather 




Figure 3: Left: The marginalized posterior probabilities of wq (up- 
per panel) and e (lower panel). Right: Icr and 2(T mean-likelihood 
contours in the (e, n'o) plane, after marginalization over all other 
parameters. Red is SNe only, blue is CMB only, and black is 
CMB+SNe+BAO. 



between px (that, for wx — -1, plays basically the role 
of no) and pz ~ H^, which has a different EOS. 

We have performed a detailed comparison of our 
model with current observations of CMB, SNe la and 
BAO using modified versions of CAMB and CosmoMC 
[37, 38], treating perturbations in the DE by model- 
ing it as a perfect fluid (without anisotropic stress) with 
the EOS parameter woe given in (12) and a unit rest- 
frame sound speed. Full details will be reported in [20]. 
In fig. 3 we give a sample of our results. The plots 
on the left are the one-dimensional posterior probabil- 
ities marginalized over all parameters except wq or e, 
respectively, while the plot on the right shows the 2- 
dimensional posterior probability marginalized over all 
parameters except the pair (e, wo). In particular, we 
find the marginalized limits -1.25 < wo < -0.908 and 
-0.0201 < e < 0.0460 at 95% C.L., consistent with 
ACDM. The means of the marginalized posteriors are 
at <wo) - -1.07 and <e) - 0.0104 and the standard 
deviations are cr^.^ = 0.0873 and cr^ = 0.0167, respec- 
tively. Thus, our wZCDM model is consistent with cur- 
rent data, and its deviations from wCDM, expressed by 
the parameter e, are constrained at the level 0(10^^). 
Future data will be able to set more stringent limits or 
to detect a non-vanishing value of e. 
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